Quantum Theory of Spontaneous Emission from Exciton-Electron-Phonon Complexes 
in Solid: Quantum Interference and Many-Body Effect 
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A full quantum mechanical theory for the spontaneous emission from excitons simultaneously 
coupled to electronic excitations and anharmonic phonons in solid is developed. Origin of detailed 
structures, such as zero-phonon line splitting, Fano lineshape near one phonon sideband, strong 
second-order phonon Stokes line, asymmetric phonon anti-Stokes lines, and two-electron satellites 
as well as their phonon replicas recently revealed in the low-temperature photoluminescence of ZnO, 
has been identified by quantitative calculations from the theory. 
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Exciton, a bound electron-hole pair which does not 
carry electric current but performs various correlated mo- 
tion of electron and hole, is of fundamental importance 
in optical properties of solids. [2, 0, 0| Over 70 year de- 
velopment has made the exciton physics become a highly 
diversified field of science and is still being actively pur- 
sued in many directions since excitons have been found 
in all types of non-metallic solids as well as in certain 
rare earth metals. The interactions of excitons with each 
other and with other quasi-particles such as phonons 
drastically affect various aspects of the exciton related 
phenomena. Furthermore, the subject is intimately con- 
nected with many other fields of physics. For example, 
the neutral-donor-bound excitons in solids have been in- 
vestigated as one of promising building blocks for future 
quantum computer. [J] 

It has been well established that the coupling of ex- 
citons with phonons leads to the phonon sidebands in 
absorption or emission spectra which reflect frequencies 
of phonons and the coupling strength Q . On the other 
hand, interactions of excitons with electronic excitations 
can also be revealed by lineshapes and radiation frequen- 
cies in the photoemission spectra 7, 8, 9]. Theoretically, 
the influences of interactions of excitons with phononic 
and electronic excitations are usually treated separately 
by perturbation theories together with Green function 
technique 

Ull ]. I n rea l solids, however, excitons are 
often simultaneously coupled to surrounding electronic 
and phononic excitations [12j. There are still difficulties 
to theoretically handle the combined effects of interac- 
tions with both phonons and electrons. Especially, the 
quantum coherence recently discovered in spontaneous 
emissions eagerly demands a firm theoretical identifica- 
tion. 

In this Letter we develop an approach to calculate 
the spontaneous emission spectrum of excitons interact- 
ing with electrons and anharmonic phonons on the ba- 
sis of exact solutions of excitonic many-body states. By 
this method the energy distributions of different exci- 
tations, the effect of quantum interference among com- 



ponents, and the many-body effect in exciton-electron- 
phonon coupling systems, are taken into account on an 
equal footing. It is found that detailed spectral features 
can be quantitatively explained by the theory. 

We consider one quantum system containing subsys- 
tems of exciton, phonon and electron. In uncoupled 
case they are expressed respectively by H ex = eoova, 
H ph = Uu) %bo, and H e i = J2i e i c l c i- Here, a, b 0) and 
Ci are annihilation operators for exciton, phonon, and 
electron at the ith level, respectively, while eo, hujQ, and 
€i are their energies. The interactions of exciton with 
phononic and electronic excitations can be written as 



Hi = a)a[V x {bl + b ) + V 2 &S&J + V^&J&J + Mo) 



+ 



i,3 



Ji,j C i C j]- 



(1) 



That is, the emergence of an exciton alters both poten- 
tials felt by electrons and by vibrations, respectively. For 
electron states the change of potential leads to violation 
of orthogonality expressed by matrix elements g^j. For 
phonons such violation may include anharmonic effect 
and is generally described by the linear element V±, the 
diagonal and off-diagonal quadratic elements V2 and V3. 
It should be pointed out that the coupling strengths gij 
and Vi t 2,3 can be tuned by the density of excitons, or by 
the power of optical pumping. 

In the initial state after pumping an exciton, we are 
working in the Hilbert subspace with one exciton and op- 
erator eft a can be replaced by 1. Thus, the phonon and 
electron parts in the Hamiltonian are decoupled. For the 
phonon part, we use the following orthogonal transfor- 
mation 



bo = ub + vtf + 
bl = utf + vb + 



w. 



(2) 



where b and b^ are annihilation and creation operators in 
the new representation obeying the boson commutation 
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rules if u 



1. In this representation the Hamilto- 



nian is diagonalized when u, v and w satisfy equations 

V 1 



Ul 



€p T ff . / 65 



2e r 



(3) 



where e„ = Kloq + V2, e r = 



4K 2 , and Z = ± repre- 



_p ; i/wu 1 r ^ , u r y tp "X f 3 

sents two irreducible sets of renormalized phonon states. 
The phonon Hamiltonian becomes 



hujibjbi + Ei 



where 



and 



E± = - 



hu; ± = (el±Wi)/e r , 



(e p + 2V 3 ) 2 2 



2e r 



(4) 
(5) 

(6) 



The elements of transformation matrix between {bi} and 
{b } representations denoted by 



Vl,n;n' = (0 



{k) n {bir' 



o) 



y/n\n'\ 

can be calculated by the iteration relation: 



(7) 



\fnui 



Vl,n-l;ri '-1 + 



Wi 



Vn + Ivi 



-Vl,n+l;n'-l 



(8) 



"T , — -T]l,n;n' — li 

with a restriction condition 

Vi,n-,n< = for n' < n, n < 0, or n' < 0, (9) 

and an initial value n/^o determined by the normaliza- 
tion condition J2„> l"J,o ; n'| 2 = 1 - 

The electron Hamiltonian can also be diagonalized 
with an orthogonal transformation Cj = Y]j, dj', 



(10) 



where e' with 



1,2,..., arc the renormalized electron 



levels. The elements g^j' and energies e- depend on {e^} 
and on Consider a general electron level struc- 

ture including hydrogen-like donor levels and the empty 
conduction band described by 



-f£, for i < M, 

S e (i — M — 1), for i>M+l, 



(11) 



where eh is the binding energy of the hydrogen-like donor, 
S e is the energy spacing in the conduction band which is 
infinitesimal but set to be finite in numerical calculations, 



and M is the number of discrete levels. Considering that 
the coupling strength is proportional to the overlapping 
of the wave functions, we adopt the following expression 
for the matrix elements 



9i 

9i% 



(ij)~ 3/2 , for i,j < M, 

j-3/2^ f Qr i < M and j > M ^ 



(12) 



where g\ and gi are adjustable parameters describing the 
coupling strengths of a hydrogen-like state with another 
one and with a state in conduction band, respectively. 

Finally the exciton-electron-phonon many-body wave 
functions can be written as 



(13) 



where |0) denotes the vacuum and n is the number of 
phonons. This state has eigen energy 



£i,z,™ = e o + Ei + nhoji + 



(14) 



The statistical distribution Fi^n^sT) of *&i t i, n de- 
pends on the non-equilibrium pumping and relaxation 
processes. Under the condition of local equilibrium, the 
distribution function is 



7^ /m\ 1 I Ei I n 

^(T^-expf- — 



(15) 



where Z is the partition function 



Z = ^2 exp 



l,n\ 



k B T ) 



Light-exciton coupling described by Hr = 
J2k W(f\a + H.c.) with f\ being the creation op- 
erator of photon k leads to annihilation of exciton 
part in state and to creation of a photon. If the 

transition begins at t = 0, the time dependent initial 
wave function can be written as 



-iEi,!,„t-rt 



(16) 



where T describes the decaying process of state "J/i^n due 
to the annihilation of the exciton. The final state with a 
photon at mode k is 



jjin (,\ _ flfio) c i' „-i(hu k +n'hu +e t ,)t 

vn! 



|0>, (17) 



where u>k is the photon frequency. The optical matrix 
element between VJ" n (*) an< ^ ^fc't™ n> (*) ^ s 



a iM ,, n , ex ^l} n n {t)\Hi\^^ n ,{t)), (18) 

where the line over the term on the right hand side de- 
notes the averaging over time. Note that the relative 
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phases of these components are coherent and can be in- 
herited by the emitted photons, the resultant spectra can 
exhibit interfering structures. So we have 



Vi, 



n\k,i' ,n' 



i[huj k - (E iilin - n'hujQ - ej/)] - r ' 



(19) 



and the spontaneous emission intensity can be written as 



p(c,r) = ^^Fu,„(T) 



E 



S(u k - oj). 
(20) 

In Eq. (|20[) effects of phononic and electronic excita- 
tions and possible coherence between them are reflected 
as illustrated in the left panel of Fig. 1. As an appli- 
cation we calculate the emission spectrum of ZnO and 
the result is shown in the right panel of Fig. 1. By 
comparison with the measured spectrum for ZnO [12j , a 
global agreement between the theory and the experiment 
is achieved. The main peak corresponding to the contri- 
bution from the components with zero phonon and with 
the lowest electron state (n 1 = 0,i' = 1) in a group of 
low-lying many-body states {^iev,i=2,n=a} is denoted as 
A. Its first-order and second-order phonon replicas are 
C and E, contributed from the components with indices 
in' = = 1) and (n' = 2,i' = 1), respectively, with 
energy separations exactly equal to Huiq. In this group 
of many-body states, the energies of the electron part is 
close to the lowest one, i.e., e' ieV — < fc^T, so the 
statistical weights from them are not negligible. Another 
set of the zero-phonon and one-phonon Stokes lines, also 
with exact energy separation tiuj , are formed from the 
components with indices i' > 1 and n' = 0, 1 in the same 
group of many-body states, and indicated as B and D in 
Fig. 1. This set is associated with electron excitations 
from the lowest states to the excited states including the 
states in the conduction band. From Eqs. (|11I19I20P 
we can see that set B-D is below set A-C by an energy 
of the order of eb of the hydrogen-like states. On the 
high-energy side of the zero-phonon line (ZPL) there ap- 
pear some weak anti-Stokes lines G and H originating 
from many-body states with higher energies. In previ- 
ous theories they are usually attributed to the radiation 
accompanied by absorption of phonons. 

In the following we discuss several important charac- 
teristics in the spectrum. If the coupling to electrons and 
the anharmonicity of phonons are set to zero, the the- 
ory reduces to the well-known Huang-Rhys theory [l3j 
as shown in Fig. 2(a), giving an important support to 
the present treatment of interactions. When the exciton- 
electron couplings are not zero, new structures appear in 
the spectrum. In Fig. 2(b) we plot the spectra of two 
cases: (i) g\ ^ and g 2 — (dotted line); (ii) gi ^= 
and gi ^ (dashed line). In case (i), the exciton is only 
bound to localized electron excitations associated with 
impurity levels. The most pronounced spectral feature is 
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FIG. 1: (Color online) Left panel: Illustration of transitions 
from a many-body state with an exciton to various final elec- 
tron and phonon states without exciton. Right panel: Cal- 
culated spontaneous emission spectrum of ZnO taking into 
account both the phonon and electron excitations. The pa- 
rameters are: e = 3.3676eV, Twj = 0.072eV, Vi = 0.012eV, 
V 2 = O.OUeV, Vz = 0.02eV, gi = -0.022eV, g 2 = 0.00036eV, 
e b = 0.038eV, S e = 0.0005eV, M = 6, Y = O.OOOleV, and 
the number of included conduction band states is 300. The 
equivalent temperature for the local equilibrium distribution 
is ksT = 27.4meV, which may be higher than the environ- 
ment temperature due to the pumping. 



the appearance of the so-called two-electron satellites on 
the lower side of the ZPL as seen from the peaks of the 
dotted curve in energy range of 3.32 — 3.34eV. The posi- 
tions of these satellite peaks depend on both the donor 
energy levels and the coupling strength with the exciton. 
Another effect of g\ is to cause the splitting of phonon 
lines including the ZPL in the Huang-Rhys spectrum. 
This splitting is often observed in fine experimental spec- 
tra as shown later. In case (ii) a coupling to the con- 
tinuous excitations associated with the conduction band 
is further included. Such a coupling produces a contin- 
uum background in the spectrum as can be seen from 
the dashed line in Fig. 2(b). It is more interesting that 
a distinct Fano dip is created between the one-phonon 
line and the above-mentioned two-electron satellites. In 
Fig. 2(c) we show the effect of the anharmonicity in the 
exciton-phonon coupling. One most direct effect of the 
anharmonic coupling is the enhancement of the 2-phonon 
line. Another effect is the new structures appearing on 
the higher energy side of the ZPL, causing the breaking 
of the mirror symmetry of the phonon Stokes and anti- 
Stokcs lines about the ZPL. Note that the anharmonic 
coupling also gives rise to a shift of the whole spectrum 
due to the change of eigenenergies of many-body states. 
In Fig. 2(d) the temperature dependence of the Fano dip 
is depicted. It is clear that increasing temperature can 
destroy the quantum coherence, resulting in the suppres- 
sion of the Fano lineshape. This trend has already been 
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observed in experiment 0. 
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FIG. 2: (Color online) Calculated emission spectra for differ- 
ent coupling cases. Remaining parameters are the same as 
those in Fig. 1. (a) Vi = Vs = gi = gi = 0; (b) Different 
exciton-electron couplings in the case of V2 = V3 = 0; (c) 
Different anharmonicity parameters V2 and V3; (d) Evolution 
of the Fano dip by increasing the temperature. 

In Fig. 3 we present a direct comparison between fine 
experimental spectra and theoretical curves in two energy 
regions of the ZPL (left panel) and two-electron satellites 
(right panel). In both the regions a general agreement 
between theory and experiment is obtained. Especially, 
most fine structures in the experimental curves can be 
reproduced in the theoretical results. This leads to a new 
understanding of the fine structures in the experimental 
spectra. 
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FIG. 3: (Color online) Comparison between fine experimental 
ZnO spectra (red solid lines) and theoretical curves (black 
dotted lines) in two energy regions of the ZPL (left panel) 
and two-electron satellites (right panel). 



In summary, a microscopic theory for the sponta- 
neous emission of excitons simultaneously coupled to 
both phonon and electron excitations is developed on the 
basis of direct solutions of the exciton-electron-phonon 
composite states. The calculated spectrum using the 
theory is quantitatively consistent with the experimental 
spectrum of ZnO. Most fine structures in the experimen- 
tal spectrum have been identified. 
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